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We present an alternative functional renormalization group (fRG) approach to the single-impurity 
Anderson model at finite temperatures. Starting with the exact self-energy and interaction vertex 
of a small system ('core') containing a correlated site, we switch on the hybridization with a non- 
interacting bath in the fRG-flow and calculate spectra of the correlated site. Different truncations 
of the RG-flow-equations and choices of the core are compared and discussed. Furthermore we 
calculate the linear conductance and the magnetic susceptibility as functions of temperature and in- 
teraction strength. The signatures of Kondo physics arising in the flow are compared with numerical 
renormalization group results. 



I. INTRODUCTION 

The single-impurity Anderson model (SIAM) is a min- 
imal model to describe the interplay of charge and spin- 
fluctuations of an interacting impurity in a metallic envi- 
ronment, including the Kondo effect. Through decades of 
theoretical research since its first proposal^ in the 1960s 
it has been thoroughly investigated. There are exact so- 
lutions from the Bethe-Ansatz-technique^ 3 - and an ac- 
curate method to describe static and dynamic properties 
is Wilson's numerical renormalization group (NRG).— ^ 

In the last twenty years renewed interest arose in 
analyzing the SIAM. One reason is the fabrication of 
nanoscale devices, in which quantum dots are coupled 
to metallic leads. In certain cases they can be described 
by the SIAM and the Kondo effect was observed.— An- 
other reason is the development of the dynamical mean 
field theory (DMFT).^ In the latter a lattice model is 
mapped to an impurity model coupled to a dynamical 
Weiss-field bath that has to be determined selfconsis- 
tently. Therefore the theoretical challenge remains to 
develop versatile and numerically inexpensive methods 
that can describe a large class of impurity models ap- 
pearing in this context. The single channel SIAM can 
serve as a benchmark to test those methods. One ap- 
proach to tackle the impurity problems that has been de- 
veloped is the functional renormalization group (fRG). 9 
Even though the fundamental equation of this framework 
is exact most methods based on the fRG are perturba- 
tive. Hence, so far it has been difficult to accurately 
resolve the non-perturbative Kondo physics. However, 
the transparency and flexibility of the fRG can lead to 
useful applications in more complex contexts, where, for 
instance, the NRG is difficult to apply. fRG-approaches 
to the Anderson impurity model come already in some 
variety, e.g. there are variants based on a frequency 
cutoff, 10-12 on Hubbard- Stratonovich-fields representing 
spin-fluctuation a 13 i 14 and on a flowing level broadening^ 
Also, non-equilibrium situations are subject to current 



research ji^— 

Here we introduce and test another fRG approach to 
the SIAM. As opposed to the previous approaches, our 
method starts with the exact solution of a small system of 
a few sites, which is termed the 'core'. The fRG flow then 
couples the core adiabatically to a bath of non- interacting 
fermions, in a 'hybridization flow'. The main motivation 
for this approach is the following. The usual hierarchy 
of fermionic fRG equations for the fermionic vertex func- 
tions has to be truncated by neglecting the higher-order 
vertex functions, typically after the four-point vertex. In 
the usual context without bare higher-order interactions 
and in standard perturbation theory, these higher-order 
vertices would come in higher orders in the bare inter- 
actions. Hence, the expectation is that the truncation 
can only be good at weaker interactions. For a normal 
many-fermion system with a full Fermi surface, in the 
beginning of the fRG flow, the higher-order terms are 
suppressed by these higher orders of the bare interac- 
tions, while at low scales, near the Fermi surface addi- 
tional phase space arguments may limit their impact. For 
strong initial interaction no argument can be given that 
the impact of these neglected vertex functions is negligi- 
ble. Another expectation is however that these higher- 
order terms are mainly determined by local physics and 
by degrees of freedom over a larger energy range in terms 
of the free Hamiltonian. Therefore one may hope to ar- 
rive at a satisfactory description also for stronger inter- 
action by incorporating the higher-order vertices of only 
a small system and by neglecting their change when the 
low-energy physics is altered during an fRG flow. Hence, 
in the present approach, we use the exact four-point ver- 
tex and self-energy of a small system as starting point 
for the hybridization flow. These quantities have built 
in the effect of all orders in the interaction at least for 
this small system. Now, performing the truncated fRG 
flow, the hybridization-induced change of the back-effect 
of the higher-order interactions on the four-point and ul- 
timately on the self-energy will be missing, but this may 
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still be better than ignoring the higher-order physics com- 
pletely. 

Note that this strategy which we are testing here for 
an impurity problem, could also be extended to a lattice 
problem. One can imagine using the small-cluster self- 
energy and four-point-vertices as initial condition for a 
flow in the band width or hopping amplitude of a lattice 
dispersion. Similar strategies have already been pursued 
for bosonic problems . 19 i 20 In this context, the present 
study can be seen as first step in the exploration of such a 
procedure for fermions, with the benefit that in impurity 
models quantitative benchmarking is possible. 

The application of RG flow equations usually requires 
a controlled starting point in the parameter space of the 
theory where the vertices are well known. Then one can 
follow the flow toward a nontrivial physical point of the 
theory. In our case the flow takes place in the effective 
theory of the first bath site next to the impurity or cor- 
related core system. Initially, the bath is decoupled, and 
hence the bath site is non-interacting, providing a well- 
defined starting point with finite density of states at low 
energies. Then, in the RG flow, the coupling to the cor- 
related core is switched on and increased to the desired 
value. Thereby the bath theory becomes interacting and 
the spectrum of the bath sites is modified. Employing 
exact relations of the first bath-site self-energy to the 
self energy of the correlated core, we can then deduce 
the spectrum of the correlated core as well and study the 
signatures of Kondo physics. 

This paper is organized as follows. In Section [TT1 
we describe the single-impurity Anderson model and its 
Green's functions on the correlated site and on neighbor- 
ing sites. In Sec. IIII1 we derive the effective bath theory 
on the first bath site and give relations between the bath 
self-energies and the self-energy of the correlated site. In 
Sec. IIV1 we describe the fRG scheme in the effective bath 
theory. Sec. |V]is the devoted to numerical results. We 
conclude with a discussion in Sec. IVII 



II. THE SINGLE IMPURITY ANDERSON 
MODEL 

A. Hamiltonian 

The Hamiltonian of the single channel SIAM consists 
of three parts 



onsite-energy is given by 



H = Hdot + Hbath + Hbath-dot ■ 



(1) 



Hdot describes the interacting electron level and is given 
by 

Hdot = E ( e d,° - aO 4<k + u 4°mK • ( 2 ) 



The operators and d a create and annihilate electrons 
on the dot-level with spin-component a = ±1. The 



U 



td.o - M = ~~2 +V g + Ba . 



(3) 



including a magnetic field term B = gp^H with Bohr 
magneton and a gate- voltage energy V g . The term 
— ¥ is choosen such that V g = corresponds to the 
particle-hole-symmetric point. 

Our bath consists of two semi-infinite tight-binding- 
chains with hopping-parameter t 

oo 

H bath = -t ]r J2 E ('v.a-i— + H - c ) 



s=L,R (7 j=l 
oo 



/< E EEO* 

s=L,R a j = l 



(4) 



The operators bj a s and bj, a , s create and annihilate elec- 
trons on the site j of the left (s = L) or the right (s = R) 
bath with spin-component er. Different choices of the 
bath could be easily incorporated into our formalism. 
The coupling between the bath and the interacting dot- 
site is given by 



Hbath-dot = -v E E { b i,<y^ d o- + H.< 

s=L,R a 

We can do a unitary transformation 



bj^a^even 
bj.a^odd 



±( 1 1 U h*,L 



(5) 



(6) 



such that only the even combination coupled to the dot 
site, since the left and right part of the chain possess the 
same chemical potential, 



Hi 



bath — dot 



-eOl 



, , en d a + H.c 



(7) 



where v = \^2v. Hbath remains formally unchanged. In 
the following we ignore the decoupled odd bath and skip 
the index 'even' on the remaining even bath. Further- 
more we set [i = 0. 



B. Green's function of the SIAM 

The free Green's function on the dot-site is given by 

G° (i(j n , d, d) = [iu n - td.o- - A (zwjp 1 (8) 

where A (iu) n ) is the hybridization function, which is 
given by 

A (iu} n ) = v 2 g b (iu} n ,bi,h) (9) 
with g b (iu) n ) = g b (iu) n , 

g b (ioj n ) = - isgn(w„)^4i 2 - (iuj n ) 2 ^j . (10) 



3 



Details for the derivation are given in appendix [A"l 

The retarded Green's function on the first bath site is 
given by 

g b (w + i0+) = ( ~ iV^ 2 - (2* - |w|) 

+ w - v 7 ^ 2 - 4t 2 9 (|w| - 2*) sgn(w)) . 



The density of states on the first bath site is then semi- 
elliptic 



p b (ui) = Im [g b (lj + i0 + )] 



M), (ii) 



with band width W — At. In most studies of the SIAM 
in the literature one considers a constant DOS and the 
wide-band limit, i.e., W is much larger than all other 
scales for the problem. Then the physics for the sym- 
metric model SIAM mostly depends on the ratio of in- 
teraction scale U and the hybridization scale A. Here 
we keep the w-dependence of the hybridization function. 
We define the quantity A = Ttv 2 pb (0) = ! r-. We choose 
for simplicity v — t, so that A = t. This means that we 
do not have two independent parameters for bandwidth 
and hybridization as one usually does for studies of the 
SIAM and the finite bandwidth actually enters the prob- 
lem. Therefore, our results differ quantitatively from the 
the wide-band limit, which is common in the literature. 
In some sense it is more similar to the first iteration of 
a DMFT calculation with a semi-elliptic DOS. We take 
Ao = t = 1 as reference energy scale in the following. We 
would like to emphasize again at this point that our for- 
malism can also deal with more general bath functions, 
as long as they can be mapped onto a linear chain with 
certain onsite and hopping parameters. 

Because the bath is noninteracting, the self-energy 
is local on the dot site. By the Dyson equation the 
full Green's function is related to Q a (ico n ) and the full 
Green's function on the dot reads 



III. EFFECTIVE THEORY FOR THE BATH 
A. Integrating out the 'core' 

We now separate the system into two parts as illus- 
trated in Fig. [TJ One part (called 'core' in the following) 
contains the correlated site and the first L bath-sites of 
the noninteracting tight-binding chain (L = 0,1,2,3). 
The other part (called 'bath') contains all bath-sites of 
the tight-binding-chain with index i > L + 1. In the fol- 
lowing we integrate out the 'core' in a functional integral 
representation of our model leading to an effective theory 
for the bath.— 




t)L+1 bi_+2 

bath 



FIG. 1. The semi-infinite tight-binding-chain is separated into 
two parts: The 'core' includes the correlated site and L bath- 
sites. The 'bath' consists of the remaining bath-sites. The two 
parts are coupled by a hopping-term which is proportional to 
the parameter A. 



Our model is described by the grandcanonical partition 
function 

Z = Jv[b,b] exp [-S (b, b)] (13) 

with the action 

S (b,b) = S core (b ,b ,bi,h, ...,b L ,b L ) 

+ Sbath f>L+l, l>L+2, &L+2, ■•■) 

(bL,bL,b L +i,b L +i) ■ (14) 



S A v 

^coupling 



To make the notation more compact, we denote the dot- 
site d as the 0th bath-site bo. We define the scalar prod- 
uct 



Go- (iu) n , d, d) = [icj n - e d , c - S djCT (iuj n ) - A {iu) n ))~ 



(12) 



Score 

and Sbath are given by 



(15) 
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S C ore{bo,bo,bi,bi,...,b L ,b L ) = (bo,( d T + £d,a)bo) +U / dr n-j- (r) nj. (r) (16) 

Jo 

L L-l 

-w [(60,61) +^.c] +^(6^6,-) -t ^ [(M*+0 + /f.c.] , 
3=1 3=1 

OO OO 

Sfcaifr (^L+l, ^L+l, 6 i+2 ,...) = X! E [(63,63+1) (17) 

j=L + l 3=i+l 

where we introduced n CT (r) = 6o,o-(t)6o !(7 (t). The coupling between the core and the bath is described by 

ScoupUng (pL , H , h+1 , b L+1 ) = -At [(b L ,b L + 1 ) + H.C.] . (18) 

I 



We introduced the flow parameter A. The original model in Eq. (fT8jl . 
((T|) corresponds to A = 1 and for A = core and bath 

are decoupled. In the case L = 0, S core is just given by The generating functional for the connected Green's 

the first line of Eq. (|T6)) and one has to replace At by Av functions of the core-problem is given by 



Wcore (J, J) = In 





$core c) 


J V [c, c] exp 


■^core J 





(19) 



with the core-fields c = (&o, 6i, 62, 6^). Z core is the partition function of the core-problem given by 

Z core = y 2? [c, c] exp [-.Score (c, c)] . (20) 



yV core can be expanded in the fields 
" (-1) 



,/3 



(—1)™ f p f p V V 

yV c .ore (J,J) = E — \T~ H / dT i- / dT " / dr i- / <*7^aSfe ) (ii,Ti;...;i„,T n |ii,7l;...;i; ) 
„ =0 n! i, f _ Jo Jo Jo Jo 



such that 



Score (H, T 1i •••> T n|*l7 T l ) •••) T n) — 



X 4 (n) . . . J in (r„) J i; «) ... J i; ft ) , (21) 
<5 2 "W COI . e (J, J) 



<yj il (Tl)...*J iB (T n )*Ji;«)...*J ii (7l) 

I 



(22) 



J=J=0 



With the definition of the (L + l)-component field yielding an effective action for the bath 
X = t(0, 0, &l+i) we can rewrite the action (fT4|) as 



S = Score (C, C) - ^ (Cj, A*.) - H .C. 



i=0 



+ Sbath (6i+l; 6i+l, &£+2j •■•) 

Now we can formally integrate out the c-fields, 

L 



(23) 



J T> [c, c] exp 



-S core (C, c) + E ( C *> A X») + H.C. 



exp [yV core (Ax, A%)] 
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Z = J V[c,c, b L+1 ,b L+ i, ...] exp 



'Score 



i=0 

= J V \pL+l,b L +l,b L +2,bL+2,.~] exp -5^, ^L+2, &L+2, ■• 

When we expand W cor e (Ax, A%) in the x, X~ ne lds the effective action has the following form 
^bath (bL+i,b L +i,b L+2 , b L+2 , ...) 

= S ba th (bL+l,b L +l,b L+2 , b L+2 , ...) - Wcore (A*, A%) 

~ (-l)"A 2n 

n=0 



+2, ■ 







^ / dn... / g?t„ / dr{ 



/3 



/3 



dr' 



x ^"i (n, n; i n , T n ri; , i' n , t£) Xh (n)-Xi„ ( T « )x< K) ■ ••Xi; ( r i ) ■ 



In the following we neglect the term with n = 0, which 
does not contain any fields. Furthermore, we truncate 
the sum over n after n = 2. This means we consider only 
the first and second order of the expansion and neglect 



all correlation- functions Gcore. ■ When we transform the 
action to Matsubara frequencies we get 



Sbath ( b L+l, b L +l,b L +2, bL+2, •■•) 

A 2 

Sbath (bL+l,b L +l,b L+ 2,b L+ 2, ...) + -g22 Xii( i<J )G%£e ( iu} , h, i[) Xi[ 
A 4 

X! Xii(«^i)x<2( iw 2)^S. ( r 2 2 (iwi,ii;iw2,i2|^^ii;«W2,^)xii(iw0Xii(w2)<ya> 1+W2) wj+u 



4£ 



3 



^ ^ bi+i,<ri (iwi)6i+i,o- 2 {iu2)Gcor} ( iuj i, b L , cry, iw 2 , &z,, o^]^, 6z, cr[;iu; 2 , b L ,a 2 ) 



Sbath (pL+1, b L+1 ,b L+2 , b L+2 , ...) H — ^ X! & i+l.°'( ia; )^cire ) ,a b L , b L ) &L + l,cr( 



ID O 

jAty 

4/3 3 

x &Z,+l,^(^l)&i+l,(r^(^2)^ 1 +W2,w' 1 -l-w^(7 1 +o- 2 ,ffi+o 



In the effective-bath theory there is a local interaction on bath site L + l. The other bath sites (L + 2,L + 3,...) remain 
noninteracting and can be integrated out. This leads to the local effective action 

S blth (b~L+i,b L +i) = -^^2^2b L+ha (iuj) (ioj- {Mf G^rla (iuj,b L ,b L ) - t 2 g b &l+i,<m" 



(At) 4 



^2 X! bL+i^ii^bL+i^ii^Gcore , b L , <t x ; iu 2 , b L , a 2 \iJ 1 ,b L , a[;iuj 2 , b L ,a' 2 ) 



3 

'"^'"i "i'"2 



4/3 



(24) 



c ( 1 ) 

The local correlation functions of the core, Gcore and 
Gcore can be calculated from the Lehmann representa- 
tion, which is given in appendix [Cl Note that for A = 0, 
the bath theory is noninteracting. The exact solution of 
this serves as a initial condition for the fRG flow in A. 



B. Relation to the Dot self-energy 

In the effective theory ([24|) the bath site L + 1 is now 
interacting with a frequency dependent term, while in 
the original theory (Tl4"f it was noninteracting. The self- 
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energy and all higher irreducible vertex-functions are lo- 
cal on the dot-site by construction . Nevertheless the lo- 
cal Green's function of the coupled problem on bathsite 
L + 1 is nontrivial and depends on the dot self-energy, 
as can be seen in (|A8|) or (| A9|) . This Green's function 
can also be derived in the setup of the effective the- 
ory and one can use the identity Q a (iu n , &l+i) = 
Q% is (iu> n ,b L+ i,b L+1 ) with 

gff (iu n ,b L+1 ,b L+1 ) = (25) 
iu - (At) 2 g c c £% (ioj, b L , b L ) - t 2 g b (iu) - S 6>(T (ia;)] 

to get a relation between the dot self-energy and the ef- 
fective self-energy £& on bath site L + l. These relations 



depend on L and for L — 0, 1,2,3 we get the relations 
given in Eq. (|B~l"j) - (|B"4| in appendix [Bj 

IV. FUNCTIONAL RG FLOW EQUATIONS 

The effective action in Eq. (|2"4")l reduces to a nonin- 
teracting model for A = 0, because the interaction term 
is proportional to A 4 . This represents a simple starting 
point for a fRG flow in A. In order to use the fRG for- 
malism for one-particle irreducible (1PI) vertices ^ 22 ' 23 , 
the flow parameter A should only occur in the quadratic 
part of the action. This can be achieved for any A ^ by 
rescaling the fields b^+i — > &l+i/A and b~L+i — > &l+i/A. 
This leads to the quadratic part of the effective action, 



S bL f h { b L+i,b L+1 ) = ~^^2^2b L+ha (iu}) -t 2 g c c ^l a {iu:,b L ,b L ) - ^gbibub^iu) 



IU) 



and quartic part which does not depend on A anymore. 
Note that the rescaling changes correlation functions of 
different order in the fields differently, but in the end we 
will study the case A = 1. The one-particle-irreducible- 
(lPI)-vertex-functions on scale A can be calculated by an 
infinite set of exact flow-equationSi 9 ' 22 i 23 If we neglect the 
flow of the three-particle- vertex and of all higher vertex 
functions we get a closed set of equations for the self- 



energy £ A and the two-particle vertex-function T A 



b ■ 



dA 



££(fc';fc) = -Tr[S A r£(fc',.;fc,.)] 



(26) 



^-r A (k[ , k' 2 ; h , k 2 ) - Tr [S A T A (k[ , k' 2 , .; k x , k 2 , .)] 
-Tr [5 A r^(., .; fa, k 2 )[g A ] T r A (k[, k' 2 ; ., .)] 

-Tr[S A Tt(.k[,.;k 1 ,.)g A Tt(k 2 ,.;k2,.)] 

-[k[ o k 2 ] - [ki O k 2 ] + [k[ o k 2 , ki O fc 2 ], (27) 

in which g A is the full propagator and S A is the so called 
single-scale propagator defined by 



(28) 



The k\ denote one-particle-quantum-numbers (in our 
case Matsubara-freqency, spin and site-index) and the 
trace is defined with respect to these quantum numbers. 
In the case B = 0, the flow-equations (|26|) and (|2T[) can be 
further simplified by using the spin-rotation-invariance of 
the effective action (|26p . which is described in appendix 

El 

By integrating the flow-equations from A = to A = 1 
we can derive the self-energy of the effective bath theory. 
Using the relations (|B1[) - (|B4[) we obtain the dot self- 
energy E d (iu). 



In the simplest approximation one neglects the flow 
of the two-particle-vertex and integrates equation (|26"j) 



with T A = r A_0 (in the following called 'approximation 
1'), where r A=0 is given by Eqs. (|DlJ and ([D9]l . If 
one integrates the full set of Eqs. (|26l) and (|2T|) (called 
'approximation 2') the numerical effort scales with the 
third power of the number of Matsubara-frequencies. 

Motivated by the fullfillment of Ward-identities in the 
RG-flow, the following replacement in the flow-equation 
for the vertex-function was proposed^ 



dg A 

dA 



s 1 ' 



G ~~d~A 



(29) 



This replacement is used in all following calculations. 

Instead of doing the fRG-flow in the effective bath- 
theory P5| one can also derive flow-equations for the dot 
self-energy (IB1|) - (|B4I) and the two-particle-(lPI)-vertex 
on the dot-site with the core-(lPI)-vertex-functions as 
initial condition. For L > these flow-equations have 
a more complicated structure than in our case and the 
calculations are easier in the setup of the effective bath- 
theory. For L = and approximation 1 we compared 
both schemes. It turned out that neglecting the flow of 
the 2-particle- vertex of the effective bath-theory leads to 
better results than doing an analogous approximation for 
the 2-particle- vertex on the dot. 



V. NUMERICAL RESULTS 

In the following we present our results for different 
sizes of the core (L = 0, 1 and L = 3). The data shown 
here was produced by integrating Eqs. (|26l I27[) numer- 
ically from A = to A = 1, using typically 100 to 200 
Matsubara frequencies at temperatures varying between 
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P = 20/A and /3 = 50/A . We set A = 1 giving the 
energy scale and in most cases consider particle-hole sym- 
metry, £d = —U/2. The results are compared with other 
fRG approaches and benchmark NRG calculations. The 
latter are carried out with the same semi-elliptic density 
of states. We focus on quantities that can directly be cal- 
culated from the data on the imaginary frequency-axis. 
One of these quantities is the effective mass m* given by 



1 = 1- 



(flmEd (iu>) 



du) 



(30) 



In the Kondo-regime the quasiparticle weight z deter- 
mines the width of the Kondo-resonance and is expected 
to scale exponentially with the interaction strength. We 
furthermore calculate the linear conductance of the dot 
G = J2 a G<j given by (in the following we set H = e 2 = 1) 



G a = -7TV< 

4*. 



dwA d ^ a (uj) p b (w) - 



driF (u>) 



driF (w) 



dui 



duA^cr (w) - 



(31) 



In the second line we used that the derivative of the 
Fermi-function is sharply peaked at low temperature at 
lu = so that the w-dependence of Pb(w) can be ne- 
glected. Of course the conductance is derived by an in- 
tegral over the real frequency axis and at first sight one 
has also to perform an analytic continuation. To circum- 
vent this, we follow an approach, proposed in Ref. l25l . 
that does not require an analytic continuation. In this 
approach G follows from the formula 



G„ 



A T Y RcM 

a>0 



dQ a (iu) a ) 
doj a 



(32) 



mine 
We 



where the imaginary frequencies iCo a and the weights 
R a are defined in Ref. 25. The frequencies iCu a differ 
from the original Matsubara-frequencies and we deter- 

dQ °dCj U)c ' f rom a Pade-approximation. 
have also done the analytic continuation 
Q (iw n , d, d) — > Q (lu + i0 + , d, d) to the real frequency axis 
using a Pade-algorithm described in Ref. [26|. As the ana- 
lytical continuation of numerical data is mathematically 
an ill-defined problem we did not obtain numerically sta- 
ble and meaningful results for all parameter sets. 



which is the atomic limit result. (iuj) diverges at 
iio = 0. In Fig. [2] (upper panel) we show the ito- 
dependenc of the self-energy at the end of the flow for 
A = 1 on the dot for U = 10 A , P = 30/ A and for the 
particle-hole symmetric case, V g — 0, computed with the 
described fRG- flow in both approximations 1 and 2. In 
the calculation we included 200 Matsubara-frequencies. 
The divergence of the self-energy at iuj = has disap- 
peared, but there is still a discontinuity, which is not 
cured by the flow. This shows the flow equations are not 
able to restore the expected local Fermi liquid properties 
of the SIAM, if we start with the atomic solution. The 
height of the unphysical discontinuity becomes however 
smaller in approximation 2 compared to approximation 
1. 
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A. The case L = 

For L = the 'core' is given by an isolated dot site. In 
this case the core-correlation-functions can be calculated 
analytically^ and one can also derive analytical results 
in the setup of the effective bath-theory^ 

The initial self-energy at A = for Ed = — U/2 is given 

by 



T, d (iuj) 



U 
1 



El 



(33) 



FIG. 2. (color online) Comparison of approximations 1 and 2 
to the Matsubara-self-energy for U = 10 Ao, P = 30/ An and 
L = 0,1. For L = one gets a discontinuity of ImE(ioj„) 
between positive and negative Matsubara-frequencies, which 
is reduced if one increases the level of approximation. For 
L = 1, approximation 1 the self-energy is continuous for small 
frequencies. As shown in the inset one obtains a small step 
between positive and negative frequencies in approximation 
2. 



In the spectral density derived from a Pade- 
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approximation to our numerical data at half filling (not 
shown) two slightly broadened atomic limit peaks at 
±U /2, but no central Kondo-resonance at small frequen- 
cies appears. Hence the L — approximation fails to 
describe the screening of the local spin-l/2-moment by 
the conduction-electrons. This screening and singlet for- 
mation should develop when A is switched on, while at 
A = the local moment is unscreened. We assume that 
this strong mismatch is the reason for the non-occurence 
of the Kondo-resonance in this approximation. 

Our results are consistent with the findings in Ref. 
I27L where a superperturbation approach to the Anderson 
model is developed. In this approach a finite local clus- 
ter containing the correlated dot-site is solved exactly. 
The correlation-functions of this cluster then serve as in- 
put for an effective theory of dual-fermion-fields. Like in 
our setup, no Kondo-resonance is found when the cluster 
contains only the correlated bath-site. 

B. The cases L — 1,2 and 3 



L = 1 consists of four delta-peaks. Two of them are 
located at £1,3 = ±\ (VU 2 + 6Av 2 + VU 2 + I6v 2 ) 
± + ^xh) ' wm ch belong to excitations from the 
ground-state ([Ml) to the one-particle-state jj-\a,e) + 
(l — jjr^j |e, cr) (v <C U) and its corresponding three- 
particle-state, which is connected by a particle-hole- 
transformation. In the limit v — > they are equal to 
the atomic ±[//2-excitations of the L=0-core. When 
we switch on the coupling to the bath in the fRG-flow, 
they evolve into hybridization broadened peaks. The 
other two peaks in the spectrum of the L=l-core lie 

at e 3 ,4 = ±i (VCT 2 + 64u 2 - y/U 2 + Wv 2 ) V % U ±Sf. 
They belong to excitations from the ground-state to the 

one-particle-state \j \a,e) + (l — jfr^j \e,a) (v <C U) 

and its corresponding three particle-state. Note, that 
due to these excitations, there is a finite spectral weight 
near zero-energy already for A — that can evolve into 
a central Kondo resonance. 



In the case L = 1 the isolated core consists of an in- 
teracting site coupled by a hopping-term v to a nonintcr- 
acting bath-site. This model is still analytically solvable. 
The ground-state at half-filling is a spin-singlet-state. For 
€4 = — U/2 and in the limit v <C U and this state is given 
by 

\S = 0) = ^ (|tl, e) + | e ,n»- (1 - ^) (it, 4.) - 14, t)) 

(34) 

with energy \(-U - VU 2 + 6iv 2 ) "w' -\ - The 
first entry in | ., .) is the correlated site, the second the ad- 
ditional uncorrelated core site, e stands for an empty site. 
Now, in contrast to the case L = 0, the local-moment on 
the dot is already in a singlet state for A = 0. This is 
a much better starting point to describe features of the 
Kondo effect. As can be seen in Fig. [2] the self-energy at 
the end of the flow for L = 1 is continuous at iuj = 0. For 
L = 2 the isolated core self-energy has a similar shape as 
in the L = 0-case and for A = 1 we get a finite step be- 
tween positive and negative Matsubara-frequencies. In 
its groundstate the core carries again a finite s = 1/2- 
moment, doublet ground state, in this case which does 
not become screened when we switch on the coupling to 
the bath in the fRG-flow. This shows once more the im- 
portance of choosing a core with spin-singlet groundstate 
for an at least qualitatively correct description of Kondo 
screening in this setup. The next larger core size with a 
singlet groundstate contains L = 3 bath sites. Numeri- 
cally the calculation of the two-particle vertex function is 
limited due to the exponential growth of the core-Hilbert- 
space. We just used approximation 1 in the L = 3 case, 
because here we only need to calculate the vertex for two 
instead of three independent frequencies. 

Let us now discuss the numerical results in more de- 
tail. The spectrum of the isolated two-site core with 



It turns out that already in the most simple approx- 
imation 1 we get a quasiparticle resonance at w = 0. 
The change of the spectrum for different values of A is 
shown in Fig. [3] for U = 6A . The peaks ei^ become 
slightly broadened, but their positions does not change 
significantly. In the end of the flow (for A = 1) their 
maxima are not located at ±{7/2 = ±3Ao, the position 
usually expected in the wide band limit with a purely 
imaginary hybridization function A(cj). However, in the 
present case where the bandwidth is less than U the hy- 
bridization function A(cj) has a finite real part which 
rcnormalizcs this position. The position of the peaks, 
±3.9Ao turns out to be comparable with what is found 
in NRG calculations. The small broadening of these high- 
energy peaks is related to the fact that they lie outside 
the bandwidth (— 2Ao,2Ao) of the bath, such that the 
width is purely due to self-energy effects. During the 
flow, already for small values of A, the low-energy peaks 
£34 become broadened and a central resonance at u) = 
with Ad(uj — 0) = ^rj- emerges. Note that at the end of 
the flow, for A = 1, there are still remnants of the peaks 
£3.4, which is interpreted as an artefact of the approxi- 
mation. The same artefacts are obtained for L = 3. 

In Fig. [3] the self-energy calculated in approximation 
1 and 2 is shown. In approximation 1 the self-energy 
is continuous for small frequencies and the derivative 

lflm ^' u ^ is negative, which leads to a reduced 

u— 0+ 

width zAo of the resonance at small frequencies. As 
shown in the inset of Fig. [5] we obtain a small step be- 
tween positive and negative Matsubara frequencies in ap- 
proximation 2. This step leads to a slight broadening of 
the central resonance, which decreases with decreasing 
temperature. Therefore it can be understood as a physi- 
cally sensible finite-temperature effect. 
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< 3 




FIG. 3. (color online) L = 1 dot-spectra for several val- 
ues of A at half filling and U = 6Aq and f3 — 50/ Ao, ob- 
tained from a Pade-approximation to our numerical data on 
the imaginary frequency axis. The atomic limit peaks become 
slightly broadened and their position changes from 4.2Ao to 
3.9Ao. At small frequencies a central resonance with height 
Ad(io = 0) = l/7rAo is emerging during the fRG flow from 
A = to A = 1. 



C. Results for the effect ve mass in comparison 

In Fig. 2] we show the effective mass for L = 1 , ap- 
proximation 1 and 2 and L = 3, approximation 1 in com- 
parison with NRG data as function of the interaction 
strength U. The NRG-data is calculated at T = for a 
semi-elliptic bath-density of states. While the qualitative 
behaviour is similar, the effective mass from the fRG- 
calculations is systematically too small compared with 
the NRG data and we can not reproduce the exponen- 
tial Kondo-scale quantitatively. For interaction-strengths 

U ~ 8-9 A the Kondo-scale T K = VFy^exp (-^) 

becomes comparable with the temperature Tk ~ jj, 
which we expect to be part of the reason for the devi- 
ations from the NRG result at large values of U. Our 
results for the effective mass fall in the range of other 
fcrmionic fRG-approaches to the Anderson-mode l 12 ' 
(cf. Fig. [5]), that are calculated for a bath in the wide- 
band-limit. A direct comparison of the data needs to 
take into account the fact that as can be seen from the 
NRG-data in Fig. [5] the effective mass for a semielliptic 
density of states with finite bandwidth is in general larger 
than for a bath in the wide-band-limit. The failure in re- 
producing the exponential Kondo-scale in the effective 
mass precisely is however common to all fRG-approaches 
to the Anderson model. Note, that all fRG-approaches 
truncate the hierarchy of flow-equations after the four- 
point level. Hence, we expect that this approximation is 
the reason for this deviation. 




U/A, 



o 



FIG. 4. Effective mass for L=l, approximation 1 and 2 
and for L=3 (/3 = 30/Ao), approximation 1 in compari- 
son with NRG-data. The NRG-data are calculated for a 
bath with semi-elliptic density of states. 



NRG large bandwidth 
NRG semielliptic DOS 
fRG L=3, Approx. 1 
fRG Jakobs 
fRG Karrasch 




U/A, 







FIG. 5. Effective mass (/3 = 30/Ao) for L=l, approxi- 
mation 2 in comparison with fRG data from Ref. [l5| and 
approximation 1 in Ref. U3 . As reference data we show 
NRG-calculations for a semielliptic density of states and 
in the wide-band- limit. 



D. Results for the conductance 

Furthermore we calculated the linear conductance G 
from Eq. f|32j) . In Fig. |6]we show G as function of the 
gate- voltage V g for several temperatures and U = 8A . 
At low temperatures /3 = 50/ A we get a plateau in the 
conductance for gatevoltages between — S- and S, which 
is due to the pinning of spectral weight at the Fermi 
energy. The plateau-value is given by the unitary limit 
2Go = 2e 2 /h. For higher temperatures the conductance 
at V g = decreases quadratically with the temperature. 

In Fig. [7] the linear conductance derived in the two 
approximation schemes 1 and 2 is shown. In approxi- 
mation 2 the linear conductance for small gate-voltages 
is reduced in comparison with approximation 1. We un- 
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derstand this again as a finite-temperature effect. In ap- 
proximation 2, the Kondo peak gets narrower, i.e. the 
effective Kondo scale comes out smaller. Hence, in this 
approximation the actual temperature is closer to 
Tk as in approximationl and the conductivity shows a 
stronger finite-temperature suppression. 




FIG. 6. (color online) Comparison of the linear conduc- 
tance for U = 8A and /3 = 20, 30, 50/ A , calculated 
in approximation 1. The value G(Vc = 0) decreases 
quadratically with increasing temperature. 



a 

--«, 

a 



approximation 1 
approximation 2 




FIG. 7. (color online) Comparison of approximation 1 
and 2 to the linear conductance for U = 8Ao, /3 — 50/ Ao. 

Fig. (]8]) shows the suppression of the gate-voltage 
at Vg = due to a finite magnetic field. As shown in 
Ref. LLIJ one can extract the Kondo-scale from this supres- 
sion within a frequency-independent fRG-scheme with 
frequency-cutoff. Therefore one defines the Kondo scale 
Tk as equal to the magnetic field B1/2 that is required 
to suppress the gate voltage G(Vg — 0) to Go = e 2 /h, 
which is one half of the unitary limit. In Fig. © we 
show B1/2 as function of U. As shown the data for 
small U can be fitted to an exponential curve of the 
form aexp (— bU/ Ao). This behaviour is expected in the 
Kondo-regime. Here we find it already for these inter- 



□ B 


= 0A Q 


□ B 


= 0.05 A 




FIG. 8. (color online) Linear conductance for U = 8An, 
P = 30/ Ao and several values of the magneticfield, cal- 
culated in approximation 1 
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FIG. 9. (color online) B 1/2 as function oiU,P = 30/A , 
approximation 1, together with an exponential fit curve. 



mediate values of U. For larger U there are systematic 
deviations from exponential behaviour. These deviations 
begin at U ~ 8 — 9A , where the Kondo scale according 
to this association becomes comparable to the tempera- 
ture, Tk ~ i = -fg 1 . From our fit we get b « 0.32, in 
good agreement with the exact value b — ir/8 ~ 0.39. 2 



E. Results for the magnetic suseptibility and 
Wilson ratio 



We also calculated the static magnetic susceptibility 
which is defined by 



d((n t )-( H )) 



dB 



(35) 



B=0 
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Here (n CT ) is the average occupation of electrons with spin 
<r, which is calculated by 



(36) 



Q>0 



with the same R a and Co a as in Eq. 1321 In Fig. [10] we 
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FIG. 11. (color online) Wilson ratio R as function of U, /3 
30/ Ao, approximation 1. in comparison with NRG-data. 



VI. CONCLUSION 
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FIG. 10. (color online) Magnetic susceptibility as function of 
U, f3 = 30/ Ao, approximation 1. in comparison with NRG- 
data. 

show the spin-susceptibility in comparison with NRG- 
data. For large values of U the spin-susceptibility is ex- 
pected to be inversely proportional to the Kondo tem- 
perature Xs ~ Therefore one expects an exponen- 
tial dependence on the interaction strength. While the 
susceptibility definitely rises with increasing U, the ex- 
ponential behaviour is not found in our fRG-approach. 
A part of this deviation might again be a thermal ef- 
fect, as for U > 9Ao the Kondo temperature falls below 
/3 _1 = Ao/30 where the calculation takes place. 

In Fig. [11] we show the Wilson ratio which is calcu- 
lated as R = X 2 + 3 X with the charge-susceptibility \c — 

lim u ^o . With the relation = A , 4 — r, 

which holds in the Anderson impurity model^^ we get 
R = A „ o7r ip . Therefore we can calculate R from our 

2m* 

data of the effective mass and the spin-susceptibility. 
R = 2 corresponds to the Kondo regime, where charge 
fiuctations are completely suppressed and the charge- 
susceptibility vanishes. As seen in Fig. [Til f° r U > 8Ao 
the fRG-data come out very close to R — 2 even though 
Xs is too small in the fRG. This points to advantageous 
cancellations of errors for this ratio in the fRG. Indeed 
m* comes out too small as well. The slight decrease of 
R for U > 9A might be again due effects of finite tem- 
perature. 



We have presented an alternative renormalization 
group approach to the single-impurity Anderson model. 
The starting point is the exact result for self-energy and 
four-point vertex of a small subsystem ('core') containing 
the correlated impurity site. Then we track the evolution 
of these two functions when the coupling to the bath is 
switched on slowly. This way, the solution of the small 
isolated cluster is implemented exactly, and the flow gen- 
erates changes of infinite order in the hybridization with 
the bath. The main approximation is the truncation of 
the flow equations after the four-point vertex. In the 
present case this means that the change of the higher- 
order vertices (six-point, eight-point, etc.) upon coupling 
to the bath is not allowed to influence the lower order ver- 
tices, i.e. the two-particle interaction and the self-energy. 
Yet, the idea that led us take this avenue was that start- 
ing with the exact two-point and four-point vertex of the 
small core contains enough strong correlation physics in 
order to give a physically reasonable results. The local 
interacting physics, such as the atomic scales are well 
represented in our approach. However, the Kondo effect 
requires to describe the subtle interplay with a continuum 
of states including many energy scales and our approach 
captures this only qualitatively, but not quantitatively. 
We expect this to be less of a problem in a self-consistent 
approach, such as DMFT, where for stronger coupling 
the Kondo resonance does not survive. 

Indeed, the numerical results for clusters with an odd 
number of auxiliary sites (L — 1 and L = 3) show that 
the flow equations produce qualitatively correct results, 
whereas for even numbers the Fermi liquid behavior is not 
recovered. The dependence of the width of the Kondo 
peak with increasing U is quantitatively different from 
NRG results, i.e. the correct exponential dependence 
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is not reproduced. For larger U the Kondo scale gets FOR 732 and FOR 912. JB acknowledges financial sup- 
smaller than the nonzero temperatures for which the fRG port from the DFG through BA 4371/1-1. 
scheme is feasible. So no clear statement can be made 
regarding the large- U behavior. However, in the inter- 
esting intermediate coupling regime, the deviations may 
be tolerable. In this sense embedding this new impurity 
solver in a different contexts to describe itinerant and 
strongly coupled physics qualitatively correctly (see dis- 
cussion below) seems a viable possibility. 

In comparison with other functional RG techniques, 
e.g. those that are perturbative in U, our data end up 
in the same range, as shown in Fig. [SJ As the ground 
state for weak and strong coupling remains the same, 
it is not entirely surprising that approaches starting at 
the opposite ends lead to qualitatively similar answers. 
The quantitative agreement could however be interpreted 
further, as a measure of what is missed by the truncation 
after the four-point vertex that is common to both lines 
of approach. 

While the use of the fRG method as impurity solver 
is for most aspects not superior to the established tech- 
niques, we hope that generalizations for correlated lattice 
systems (i.e. with more than one correlated site, like the 
two-dimensional Hubbard model) will be feasible. In this 
case, both self-energy and interaction vertex will become 
increasingly non-local during the flow, and certainly, suit- 
able approximations have to be found in order to keep the 
amount of information manageable. For example, small 
correlated cluster cores can be coupled together during 
the flow via switching on the hopping amplitude between 
the clusters from to the original value. The solution of 
the core will then provide the spectral weight transfers 
on the energy scale U and the accompanying reduction 
of the spectral weight near the Fermi level. Together 
with the core interaction vertex this spectrum will serve 
as effective action of a strongly correlated Fermi liquid, 
which then can undergo long-range ordering transition 
when the cores are coupled together. Note that in ex- 
tension of earlier ideas in the vein of cluster perturbation 



lows one to determine the non-local hybridization effects 
on the interaction, which has direct consequences on the 
character and scale of low-temperature instabilities such 
as unconventional superconductivity. This way we hope 
to extend the successful functional RG instability anal- 
ysis for weakly correlated fermions to the more strongly 
correlated regime. The high-energy physics of a strongly 
interacting Hubbard-like system is certainly more local 
than the low-energy physics of collective ordering. Hence, 
the break-up into small cores and subsequent coupling to- 
gether also closely follows the physical intuition of first 
solving the problem with the largest energy scale before 
the low-energy end is considered. 

We acknowledge useful discussions with David Joerg, 
Manfred Salmhofer, Sabine Andergassen, Severin Jakobs, 
Christoph Karrasch, Volker Meden, Andrej Katanin, 
Walter Metzner, David Rosen and Walter Hofstetter. 
This project was supported by the DFG research units 
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Appendix A: Expressions for the Green's functions 



identity 



The inverse free Green's function Q (iu) n ) = 
[G° (iu> n )} on the imaginary frequency axis is given 

by Q fa n ) — WnX — Hq, where Hq is the noninteracting 
part of the Hamiltonian ([1]). Written as a matrix Q fa n ) 
it is given by 



fan) = 



' Qt ( iuJ n) 








Ql fan) _ 



where 



Qa fan) 





d 


bi 


b 2 


d 


iu) n — €d,a 


V 




h 


V 


iu>n 


t 


b 2 




t 


iUn ■ ■ ■ 



(A2) 



The inverse of Q a (iu> n ) can be calculated by using the 



" A 


B 


C 


D 



(A3) 



(A-BD^C) 1 


-(A - BD~ 1 C)~ 1 BD~ 1 


-D~ r C (A - BD- l C)~ l 


(D-CA- l B)- x 



which is valid for arbitrary invertible matrices A, B, C 
and D. 

If we introduce the matrix <?{, (iui n ) by 









h 


b 2 




(Al) 




bi 




t 






9 b 1 fan) = 


b 2 


t 




t 






b 3 




t 


iuj n ■ ■ ■ 



(A4) 



the free Green's function on the dot site follows from the 
identity (|A3I) as 

Q° fa n , d, d) = fa n - e d _ a v 2 g b fa n , h,bi)) 1 . 

(A5) 

The function gt [iu) n , b\, b\) can be calculated again from 
the identity (|A3[) . this time applied to Eq. (IA4[) . as 



g b fanM,b\) = fan - t 2 9b fanMM)) ■ (A6) 

Here we used that adding or removing the first bath-site 
from a semi-infinite tight-binding-chain do not change 
the chain. Eq. (|A6|) can be solved to give the explicit 
expression for g b fa n , 6i, foi) in Eg . (flT))) . 

By the Dyson equation the full Green's function is re- 
lated to Q a fan) and the self-energy 



[Qa fan)} 1 = Qa fan) ~ S CT fa n ) 



(A7) 





d 


h 


b 2 


d 


iuj n - e dM - Sd^cr fan) 


V 




h 


V 




t 


b 2 




t 


iuj n ■ ■ ■ 



Inverting this matrix with the identity (|A3|) gives the full Green's function on the dot Q a fa n ,d,d) in Eq. ([T2 
In the same way one gets the Green's function for the first and second bath-site, 



Qa fanMM) 



iu n - e d ,a - Sd )CT fa n ) 



Qa fanMM) = ib) n - 



t 2 g h fanMM) 



t 2 g b fa n MM 



(A8) 
(A9) 



iw„-ed,<T-S d ,„(iw„) 



For the other bath-sites with site index > 2, analogous expressions can be derived. 
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Appendix B: Relation between bath and dot self-energy 



From the identities of the Green's functions in section UlI Bl we can derive the following relations for the self-energy 
for L = 0, 1,2,3, 

(Av) 2 

L = : S d)ff (iu) n ) = iui n - e d , CT ^ , (Bl) 

core.a 

v 2 

L = 1 : £ d)CT (iuj n ) = iuj n - e d , a - ^—^ , (B2) 

1(jJ — - - 

" E 6l<r (iu, n ) + (At) 2 g&rl<r(iUn,blM) 

V 2 

L = 2 : S d)ff (iuj n ) = iui n - e d ^ - p , (B3) 



S i ,_„(^„) + (At)2 g^l a (_i^ n , b 2 ,b 2 ) 

v 2 

L = 3 : Ed.o- (iwn) = «w„ - e d , a p • (B4) 

«w n - -j 



r (i^„) + (At)2e=i ( r 1 2 i „(i^„,b3,6 3 ) 



Appendix C: Solution of the 'core'-problem 



To solve the local core-problem, we diagonalize the 
'core'-Hamiltonian exactly. From this solution we de- 
rive the local correlation functions using a Lehmann- 
representation. For the one-particle-Green's function this 



representation is given by 
(iuj n ,i,j) 



1 exp [- /3E m ] + exp [-fiE n 
7. h, 



(E n — En 



x (n\ci\m) (n\cj \m) 



(CI) 



The Lehmann-representation of the two-particle-Green's 
function is derived in Ref. It is given by 



Q {2) (i[,a,iLu[;i' 2 ,a\iuj' 2 \i 1 ,cr,iuj 1 ;i2,cr' ,iuj 2 ) = \ ^ ^(t>(Ei,Ej,Ek,Ei, iuim , i^u 2 , ^n 3 ) 

i,j,k,l n 

xsgn(n) (i\Ou,\j){j\Ou 2 \k){k\0 U3 \l)(l\ci 2 , a .\€) 



(C2) 



Here the frequencies corresponding to creation and annihilation operators have the same sign. The operators 0, are 
defined by 0\ — ct, a , C 2 = ct, cr , and O3 = c^.o-. The function <f> is given by 



4>(Ei,Ej,E k , Ei, iwi,Mj 2 ,iw 3 ) 



1 



?w 3 + E k ~ E t 



i (U>2 + W3) + - ^ 



iwi + Ei — Ej i (uii + uj 2 + 0J3) + Ei — Ei 



(iuii+Ei-Ej) iuJx + Ei-Ejl iu 2 + Ej - E k \ ioJi + E t — E 



(1 — 5 u)1 -u) 2 &Ei,E k ) ■ 



e -m _ e -f3E k 



i (uji + ui 2 ) + Ei - E k 
From (|C2p one gets the connected two-particle-Green's function from the relation 



(C3) 



gc,(2) 2 /| 1; 2 ) = g(2) ^ 2 /| 1; 2 ) _ ^ ^/|^ g(l) ( 2 '| 2 ) + £ (l'|2) 0^ (2'|1) 
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Appendix D: fRG-flow-equations 

Because the action (f26|) is spin-rotation-invariant, the fRG-flow-equations (|26|) and (|2T[) can be further simplified. 
Using the spin-conservation, the two-particle vertex is given by 

T b (iuj[ , a[ ; iw' 2 , a' 2 \ iojt , oi ; iw 2 , cr 2 ) = VJ, A (iw^ , iw 2 \iui , iw 2 )^ CTl , CT j 5 CT2 ,^ - V^ A (iwj , iu' 2 \ iut , iw 2 )(5 (Jl y 2 S a2 ^ . 

From the antisymmetry of r A (1', 2'|1, 2) under the permutations 1' O 2' and 1 O 2 it follows that the functions Vj, A 
and VJ, obey the relation 

U fc A (iwijiwjliwi,^) = Vb A (iw2,iwi|iwi,ia>2) = ^ A ("^, iw 2 |iw 2 , iwi) . 
Using this parametrization we get the flow-equations 

-^-£ A (iw) = -i ^ S A (<w') (2U b A (iw, iw'liw, io/) - U b A (zw, io/lity, iw)) (Dl) 



(iwi,tt4|ia;i,W2) = $ A p iu' 2 \iui, iu> 2 ) + $ A p/t (iu>£, w 2 |iwi,iw 2 ) + $ A p/l (iu4,iw 2 |wi,m; 2 ) (D2) 



with 



^pp (iwijiw^iwijiwa) = L (iw 3 , iw 4 ) U A (icj 3 , ii^lwi, «w 2 ) U A (iwi, kj 2 |iw 3 , iw 4 ) (D3) 



2^3 ,ICU4 



$ A p/l (iw^iwaliwiji^) = ^ L(iw 3 ,iw 4 ) ^2U A (iwi, iw 3 |iwi, mj 4 ) U b A (iw 2 , iw 4 |iu; 2 , iw 3 ) 



-U b A (iwi,w 3 |ic«;i,ia;4) Vj, A (iw 2 , iw 4 |iw 2 , iw 3 ) 

-U b A (iwi,icj 3 |iw 4 ,ia;i) U A (iw 2 , iw 4 |iw 2 , iw 3 ) ) (D4) 



$ A pft (iwi,iw 2 |icJi,ia;2) = -3 ^ L (ito 3 , iw 4 ) U A (iw 2 , «W3|w 4 , iui) V b A (iu[, «w 4 |iw 3 , iw 2 ) . (D5) 
The function i is defined as 

L (iwi,iw2) = £ A (iwi) S A (tw 2 ) + £ A (*W2) 5 A (iwi) . (D6) 
The single-scale-propagator is given by 

s A M = -2K^-e 9b ^MM)) {m) 

(iu - (tkyg c £} (iuj, b L ,b L ) - t 2 g b [iu, h, bi) - A 2 £ A (iuf 

The initial conditions for A = are 

E A =° (iw) = 0, (D8) 

U b A=0 (iwi.iwiltwi.iwa) = (twi, 6i, t; ««4 &x4 &£, t; *W2, 6 L , |) . (D9) 
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